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New Neighbours: Modelling the Growing 
Population of 7-ray Millisecond Pulsars 

C. Venter, A.K. Harding, and T.J. Johnson 


Abstract The Fermi Large Area Telescope, in collaboration with several groups 
from the radio community, have had marvellous success at uncovering new y- ray 
millisecond pulsars (MSPs). In fact, MSPs now make up a sizable fraction of the 
total number of known y - ray pulsars. The MSP population is characterized by a va- 
riety of pulse profile shapes, peak separations, and radio-to-y phase lags, with some 
members exhibiting nearly phase-aligned radio and y- ray light curves (LCs). The 
MSPs’ short spin periods underline the importance of including special relativistic 
effects in LC calculations, even for emission originating from near the stellar sur- 
face. We present results on modelling and classification of MSP LCs using standard 
pulsar model geometries. 


1 Introduction 

To date, Fermi Large Area Telescope (LAT) has uncovered more than a dozen new 
y-ray millisecond pulsars (MSPs) [17]. These MSPs were discovered by folding the 
y-ray data with the known radio period, using timing solutions provided by various 
radio telescopes. Discovery of the first eight Fermi - LAT MSPs [1] indicated that 
these MSP pulse profiles and spectral properties mimic those of the younger pulsars, 
and that there are non-zero lags between the radio and y-ray peaks. However, the 
discovery of PSR J0034— 0534 [3] exhibiting phase-aligned y-ray and radio peaks 
provided evidence for co-location of the radio and y-ray emission regions. This 
phenomena has previously only been observed for the Crab pulsar. Recent detection 
of y-ray pulsations aligned with the radio peaks [23, 4] from the famous ‘first’ MSP 
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PSRJ 1939+2 134 [6], as well as the first ‘black-widow’ MSP PSR J1959+2048 [15], 
now provide two more examples of this exceptional behaviour. Publication of a 
number of new y-ray MSP detections are awaited [22] (some resulting from follow- 
up radio observations of unpulsed Fermi-LKT sources which appear in the First 
Fermi-LAJ source catalogue), and the growing MSP population now constitute a 
significant fraction of all Fermi-LAJ pulsar detections to date [2]. 

The y-ray MSPs exhibit a variety of light curve (LC) shapes, and have been mod- 
elled using standard two-pole caustic (TPC), outer gap (OG), and pair-starved polar 
cap (PSPC) models [26], an annular gap model [11], lately also altitude-limited TPC 
/ OG models [25], and lastly using a force-free magnetic field [7]. We elaborate on 
some of these models in Section 2. (See [ 1 7] for a review of y-ray pulsar population 
modelling.) 


2 Pulsar Models 

Below is a brief description of some modelling aspects (see [26] for details). 

B -field As before, we assume the retarded vacuum dipole field [10, 13] to be the 
basic geometric structure of the MSP magnetosphere. The magnetic moment is in- 
clined by an angle a with respect to the spin axis, while £ is the observer angle. 

Caustics Due to the rapid increase of the co-rotation velocity with altitude, ap- 
proaching c as one reaches the light cylinder, it is important to include relativistic 
effects when calculating LCs. The effects of (i) B-field line curvature, (ii) time-of- 
flight delays, and (iii) aberration combine to form intense caustic emission on the 
trailing field lines, while the emission is spread out on the leading field lines [20, 14]. 

Polar Cap (PC) Model The PC model [9] assumes extraction of primary charges 
from the stellar surface, acceleration of these charges to relativistic velocities along 
B-field lines, and subsequent y-ray curvature radiation and pair creation at a very 
low altitude. A pair- formation front (PFF) is formed which screens the accelerating 
E-field (£]|) parallel to the B-field. 

Faded PC Model This is a geometric approximation [25], where we use a rising 
and falling exponential function to modulate the assumed emissivity close to the 
stellar surface, as motivated by pair cascade modelling. 

PSPC Model When £] is not large enough to accelerate primary electrons to suf- 
ficiently high energies so pair cascades will occur, it will be unscreened, and one 
expects acceleration and radiation from the full open field line region above the 
PC [18, 19]. 

SG / TPC Model When enforcing the boundary condition that £] vanishes on the 
last open B-field lines (LOFLs), the PFF altitude becomes dependent on magnetic 
co-latitude, and narrow gaps form just inside of the LOFLs [5]. The SG model [21] 
may provide the framework for the geometric TPC model [12]. 
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Fig. 1 Altitude-limited TPC ! OG LC fits for PSR J 1939+2 134 (panel [a]: '/-rays', panel [b]: 
radio) with («™ C i PC) = (65 °,80<=), {a OG !C o G) = {70 - |!W)t = 

(0.12,1.05)* LC , (R^.R^) = (0.12, 1.00)/? LC , {i?“ h n °’ 1ITC ,^ 0 ’ TPC ) = (0.93.1.20)« lc , 
(^mdi 1 ,, OG J? md,a.° G ) = ( 0 .75, 1 . lO)/?^. Fits for PSR J1 959+2048 (panel [c]r y- 
rays\ panel [d]: radio) with (a TPC ,C TPC ) - (60°, 80°), (a 0G ,C° G ) = (55°, 85°), 
(*^-^ PC ) = = (0.12,1.20)i?Lc, (/S > ’ TTC ^ 0 ' TPC ) = (0.95,1.05)^c, 

Kt°'° G - R ™ 00 ) = (0-90, 1.10)*lc. 


OG Model The Goldreich- Julian charge density [16] changes sign at the null charge 
surface (NCS), and OGs are assumed to form close to the LOFLs, above the NCS [8, 
24]. 

Altitude-limited TPC / OG Models We lastly investigated geometric TPC / OG 
models with arbitrary lower, J? m ; n , and upper, /? max , radial gap extensions [3]. 


3 Results 

First 8 MSPs We found [26] that 6 of the 8 Fermi-LKY MSPs may be modelled 
using standard TPC / OG models, while 2 are exclusively fit by PSPC models. 

PSR J0034-0534 We found altitude-limited TPC / OG LC fits [3], with (a, Q = 
(30°, 70°), (R^RLd = (0.12,0.9)fi LC , and (*Sg°,SS£>) = (0.6,0.8)* LC , with 
Rye = cP fin , and P the pulsar period. There also exist solutions [25] for the faded 
PC model for (a, Q = (1 0°, ~ 34 r ), although the first ones seem more credible. 
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PSR J1939+2134 and PSR J1959+2048 Figure 1 indicates altitude-limited TPC / 
OG LC fits with aligned radio and 7 - ray pulses [4]. 


4 Discussion and Conclusions 

Our results indicate that the high-energy (HE) radiation comes from narrow gaps in 
the outer magnetosphere, requiring the presence of robust pair formation to screen 
much of the open field line region. However, MSPs do not produce enough pairs in 
the standard vacuum dipole models to screen £|j. Standard geometrical models still 
provide suitable LC fits though, and the resulting a and £ are in reasonable agree- 
ment with fits from radio polarization measurements. In the case of phase-aligned 
7 -ray and radio pulses, we infer the existence of co-located emission regions, as well 
as high-altitude radio emission. Future work includes investigations of HE spectra 
using full radiation codes, collective emission from MSPs in globular clusters, and 
pulsar population studies. 
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